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Abstract 

In this work we present an accurate parameterization of the anti-neutrino flux pro- 
duced by the isotopes ^^^U, ^^^Pu, and ^^^Pu in nuclear reactors. We determine the 
coefficients of this parameterization, as well as their covariance matrix, by performing 
a fit to spectra inferred from experimentally measured beta spectra. Subsequently we 
show that flux shape uncertainties play only a minor role in the KamLAND experiment, 
however, we find that future reactor neutrino experiments to measure the mixing angle 
013 are sensitive to the fine details of the reactor neutrino spectra. Finally, we investi- 
gate the possibility to determine the isotopic composition in nuclear reactors through 
an anti-neutrino measurement. We find that with a 3 month exposure of a one ton 
detector the isotope fractions and the thermal reactor power can be determined at a 
few percent accuracy, which may open the possibility of an application for safeguard 
or non-proliferation objectives. 
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1 Introduction 



Experiments at nuclear reactors have a long tradition in neutrino physics. Starting from the 
experimental discovery of the neutrino at the legendary Cowan- Reines experiment [1] , many 
measurements at nuclear power plants have provided valuable information about neutrinos. 
For example, the results of the Gosgen [2], Bugey [3], Palo Verde [4], and CHOOZ [5] 
experiments have lead to stringent limits on electron anti-neutrino disappearance. Reactor 
neutrino experiments have become very prominent again due to the outstanding results of 
the KamLAND experiment [6,7]. For a review on reactor neutrino experiments, see Ref. [8]. 
Recently the possibility to determine the last unknown lepton mixing angle ^13 by a reactor 
neutrino experiment with a near and far detector is actively investigated (see Ref. [9] and 
references therein). Building on the experience gathered in oscillation experiments ideas of 
"applied neutrino physics" appeared [10-12]: A detector close to a nuclear reactor could be 
used for reactor monitoring, either for improving the reliability of operation of power reactors 
or as a method to accomplish certain safeguard requirements in the context of international 
treaties for arms control and non-proliferation of weapons of mass destruction. 

The standard detection process for reactor neutrinos is inverse beta decay: 

Ue + p + n . (1) 

The cross section a{E^) for this process is very well known [13], at an accuracy better 
than 1%. Per fission roughly 6 electron anti-neutrinos are produced (see e.g. Ref. [8]), with 
energies peaked around 1 MeV. However, for inverse beta decay only neutrinos with energies 
above the threshold of 1.8 MeV are relevant. In nuclear reactors electron anti- neutrinos in 
that energy range are produced dominantly by the beta decay of the fission products from 
the four isotopes £ = ^^^U, ^^^Pu, ^^'^U, ^'^^Pu.^ We denote the fiux from the isotope i by 
(j)i{Ey) in units of anti-neutrinos per fission and MeV. In Tabled the total number of z/g per 
fission above 1.8 MeV is given for ^^^U, ^^^Pu, ^"^^Pu, and ^^^U. 

Accurate information on the anti-neutrino fiux from ^^^U, ^'^^Pu, and ^^^Pu can be ob- 
tained by the measurement of the beta spectra from the exposure of these isotopes to thermal 
neutrons [16-18]. Subsequently these beta spectra have to be converted into anti-neutrino 
spectra, taking into account the large number of beta branches involved. These spectra 
are in excellent agreement with the direct observation of the anti-neutrino spectrum at the 
Bugey [19] and Rovno [14] reactors. The errors on these fiuxes are at the level of a few 
percent. Since modern reactor neutrino experiments aim at precisions at this level, a proper 
treatment of the fiux uncertainties becomes necessary. For ^^^U, which does only undergo 
fast neutron fission, no similar measurements exist, and one has to rely on theoretical cal- 
culations [20,21]. 

In absence of neutrino oscillations the number of positron events for a measurement time 

-'^The next to leading contributions come from the isotopes ^''"Pu and ^''^Pu and are of the order 0.1% 
or less [8]. Further sub- leading effects are the beta decay of ^^^U, ^^^Np, ^^^U (produced by radiative 
neutron capture), and corrections to the spectra from fission fragments due to neutron absorption by these 
fragments [14]. These effects are relevant for the low energy part of the anti-neutrino spectrum <2 MeV 
and will be neglected in the following. See also Ref. [15]. 
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iV; [MeV] 


235^ 
238 u 

241PU 


1.92(1 ±0.019) 201.7 ±0.6 

ooo/i 1 r\ r\c\r\\ c^rw f\ \ r\ r\ 

2.38(1 ± 0.020) 205.0 ± 0.9 
1.45(1 ±0.021) 210.0 ±0.9 
1.83(1 ±0.019) 212.4 ±1.0 



Table 1: Total number of v^, per fission above i.8 MeV (see Section |31 for details) and energy release per 
fission (reproduced from Table 2 of Ref. [5]) for the isotopes relevant in nuclear reactors. 

T in a given positron energy bin % can be calculated by 

iV. = 1^ E / ^(^-) ^*(^-) • (2) 

Here is the number of protons in the detector, L is the distance between reactor core 
and detector, and Ri{E^) is the detector response function for the bin i (including energy 
resolution and efficiencies). If the initial composition of the reactor fuel is known, the number 
of fissions per second Nf^ of each isotope i can be calculated accurately (better than 1% [8]) 
at each burn-up stage by core simulation codes. The thermal power output P of the reactor 
IS given by P = Y^i^e'Ee, where Ei is the energy release per fission for the isotope i, 
see Table ^ Since the errors on Ei are less than 0.5% we will neglect in the following the 
uncertainty induced by them. Defining the relative contribution of the element i to the total 
power^ fi, the A^^*^ can be expressed by P and Eq. (jj)) becomes 

iV, = ^ P 5^ ^ dE, a{E,) MEu) R^{E,) with /, ^ . (3) 

The aim of the present work is to consider various aspects of reactor neutrino spec- 
troscopy, with a main emphasis on issues related to the emitted anti- neutrino fiux. First, 
in Section |2l we present a phenomenological parameterization for the anti-neutrino fiuxes ipi 
based on a polynomial of order 5. We show that in many situations existing parameteriza- 
tions [22] are not accurate enough to describe the reactor neutrino spectrum at the required 
level of precision. In addition, in Section El we give a detailed consideration of the errors 
associated to the reactor anti-neutrino spectrum and provide them in a suitable form for 
implementation in data analyses. In the following we consider the implications of uncertain- 
ties of various quantities appearing in Eq. (jH} for several experimental configurations. In 
Section m we discuss the impact of errors on P, /£, and (f)i for the KamLAND experiment, 
whereas in Section El we discuss the relevance of our new parameterization of (/)^ and its 
errors for future reactor experiments to measure 6'i3. In Section IHl we consider the potential 
of an anti-neutrino detector close to a reactor: In Section 16.11 we discuss the improvement 
on the flux uncertainties from a near detector, and in Section IfT^ we investigate the possibil- 
ity of reactor monitoring by using the anti- neutrino measurement. In contrast to previous 
studies [10-12] we employ full spectral information, which allows the determination of the 
isotopic content of a reactor, i.e. the fractions /£, as well as the reactor power P without any 
external information. We conclude in Section [71 

^The power fractions fi must not be confused with the relative fission contributions ff^ = Nf^ / Nf^. 
In this case one would obtain Nf"" = ff''P/{E), where the mean energy per fission is given by {E) — ff^Eg. 
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2 A parameterization for the reactor anti-neutrino flux 

In Refs. [16, 17] anti-neutrino spectra from the fission products of ^^^U, ^^^Pu, and ^^^Pu are 
determined by converting the precisely measured associated beta spectra. In this section we 
propose a phenomenological parameterization for the reactor anti-neutrino flux, based on 
these measurements. Similar as in Ref. [22] we parameterize the spectrum of a given element 
using a polynomial: 



The coefficients au are determined by a fit to the data of Refs. [16,17]. To this aim we 
minimize the following x^-function: 



where and (jr^ = (pi^Eu ) are the values of the neutrino energy and the corresponding 
anti-neutrino flux, respectively, provided in the tables of Refs. [16, 17] for values of the 
neutrino energy By ranging from 1.5 to 9.5 MeV in steps of 0.25 MeV. Since we are fitting 
the logarithm of the flux the covariance matrix Sij contains relative errors of the For 
the diagonal elements Su we take the errors as given in the tables of Refs. [16,17] (converted 
from 90% CL to la and squared), which contain the statistical error from the beta spectrum 
measurement, a systematic error on the overall calibration, and a systematic error from the 
conversion from beta to anti-neutrino spectrum. The off-diagonal elements are obtained from 
the error on the absolute calibration, which is taken as fully correlated: Sij — cr^^'cr^^^ for 
i ^ j- The errors aj^^ are given at two calibration energies for each isotope in Refs. [16-18], 
and we interpolate linearly between these reference points. Note that this procedure assumes 
that the systematical errors from the conversion from beta to anti-neutrino spectrum are 
completely uncorrelated between different energies. 

First we have performed a fit of a polynomial of second order [K^ = 3). The resulting 
coefficients ake are given in Table |2l in Appendix and are in reasonable agreement with 
the ones obtained in Ref. [22] (some deviations appear for ^^^Pu). However, we find that 
the quality of this three parameter fit is very bad for all three isotopes (the is given in 
Figure H in the following). We conclude that at the level of precision provided by the errors 
the data cannot be described with sufficient accuracy by the polynomial of order 2. We 
have checked that a reasonable fit is obtained for all three elements only by going up to a 
polynomial of order 5, corresponding to Ki = Q parameters. The best fit coefficients ake are 
given in Table El in Appendix El and the corresponding anti-neutrino spectra are available 
in computer readable format at the web-page Ref. [23]. 

Our fit is illustrated in Figure ^ where we show the resulting spectra for the 3 and 6 
parameter fits in comparison to the data. Large differences between the 3 and 6 parameter 
fit are visible by eye only for the high energy region, where the spectra are very small and 
errors are large. However, comparing the corresponding x^-values xfa) and xfe) given in the 
figure it is obvious that the 6 parameters are necessary to obtain a reasonable goodness- 
of-fit. In the lower panels we show the residuals of the fit, i.e., for each data point i we 
plot (0£data ~ ^ilt)/'^i^ whcrc the error is obtained from the covariance matrix S hy ai = 
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Figure 1: (Color online) Illustration of the fit to the data on the anti-neutrino spectra from ^ssy j^g]^ 
239p^ [17], and ^^^Pu [17]. The thick/red curves correspond to a 6 parameter fit (polynomial of order 5), 
whereas the thin/blue curves correspond to a 3 parameter fit (polynomial of order 2). Also shown are 
the data with their la error bars and the P^r degree of freedom (number of data points minus fitted 
parameters). In the lower panels we show the residuals of the fits. Note that because of correlations the 
shown residuals do not add up to the given x^-values. The data points to the right of the dotted line in the 
^^^U-panels are excluded from the fit. 



'i^MataV^^- Note that because of correlations between the ^^^^ta these residuals do not add 
up to the total x^- Fo^' ^^^U the 3 parameter fit shows rather large residuals over the full 
energy range. Since this isotope gives the main contribution to the reactor anti-neutrino flux 
it is very important to model its neutrino spectrum correctly. Let us note that we exclude 
six data points at high neutrino energies from the fit. The change in the spectral shape 
around 8 MeV [16] cannot be fitted very well by the polynomial,^ although by accident the 
6 parameter fit gives a reasonable approximation also in this energy range. Also for ^^^Pu 
the high energy range MeV is important. In this case it is possible to obtain a good 

fit from the polynomial of order 5 even including the high energy part. 

^In Ref. [14] a term of order E"^" is introduced to model this sharp falloff. 
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3 Quantifying the anti-neutrino flux uncertainties 



In this section we discuss in detail the uncertainties on the reactor anti-neutrino fluxes. In 
TableElin AppendixElwe show the errors da^ii on the coefficients of the polynomial as well as 
their correlation matrix p^^,, as obtained from the fit to the measured beta spectra. Hence, 
the covariance matrix for the coefficients can be obtained by 

V^^, = 5aki5ak'ipik' ■ (6) 



From the table one observes that for a given element the coefficients are strongly corre- 
lated or anti-correlated, since for most elements of the correlation matrix we obtain \Pkk'\ ~ 1- 
Therefore, we perform a rotation in the space of the a^, such that the covariance matrix 
becomes diagonal. Let us for each isotope introduce new coefficients Ckt by 



^k'k Ck'e , 



(7) 



where the orthogonal matrix is defined by 

OW'iOr = dmg[i6cker] • 



(8) 



Hence, the Scke are the (uncorrelated) errors on the coefficients Cki- Using Eqs. (jH) and ((Tj) 
the anti-neutrino flux for the isotope i can be written as 



exp 



J^CkipiiE^) 



k=l 



(9) 



where pi{E^) is a polynomial of given by 



Pk 



(10) 



k'=l 



These polynomials describe the uncorrelated contributions to the error on the anti- 
neutrino flux. For example, let us consider some observable X, involving the anti-neutrino 
flux in the following way: X = J dE^ h{E^) (j)£{E^) , where h{E^) is some function of the 
neutrino energy. Then the error contribution from the coefficient Cki is given by 

dX f 

6X = 6cM dE, h{E,) ME.) ScmpUE,) . (11) 

Hence, the product 6ckipi{Ei,) is a measure for the importance of the error 6cke for any 
observable. In the upper panel of Figure |21 we show the polynomials Eq. ()10j) weighted by 
the corresponding error for ^^^U. For ^^^Pu and ^"^^Pu we obtain very similar results. One 
observes that in the relevant range of the anti-neutrino energy the flux uncertainties are at 
the level of 2%. The weighted polynomials 6ckepi{Et,) for the isotopes ^^^U, ^^^Pu, ^"^^Pu are 
available in computer readable format at the web-page Ref. [23]. Once these functions are 
known, the flux uncertainties on any observable can be included similar to Eq. pi|) . 
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Figure 2: (Color online) Uncorrelated anti-neutrino flux uncertainties for ^^^U. The upper panel shows 
the polynomials p^^{E^) given in Eq. (|f 0|l multiplied by the corresponding error 5cm- The lower panel 
shows the functions ScMPki^v) cr{E„) (l)e{E^), where a{E^) is the detection cross section. The dashed curve 
corresponds to a{E„) (p£{E^)/100. 

As first simple application let us mention how one can calculate the number of anti- 
neutrinos per fission above the threshold and its uncertainty, as given in Tabled Given 
the best fit parameters and their covariance matrix for ^^^U, ^^^Pu, ^^^Pu we readily obtain 

Jl.SMeV ^ OO'klOO'k'l ^ \OCm J 

In addition to the three isotopes ^^^U, ^^^Pu, ^^^Pu also 238u gives a contribution of a few 
percent to the reactor anti-neutrino fiux. For this isotope no measurements exist and one has 
to rely on theoretical calculations [20,21]. In the following we will always adopt for ^ss-q ^j^g 
parameterization with the second order polynomial given in Ref. [22], which we reproduce in 
the last row of Table |21 In particular, that parameterization has been used to calculate also 
the value of for ^38^ given in Tabled the error of 2% is an educated guess motivated by 
the errors obtained for the other isotopes. Since no covariance matrix of the fiux coefficients 
for 238u jg available we will always assume in the following that they are known exactly. 
Since the contribution of ^38^ to the total fiux is rather small, this assumption has very little 
impact on the conclusions drawn in this work. 
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4 The impact of anti-neutrino flux uncertainties in KamLAND 



The KamLAND [6, 7] reactor neutrino experiment is located in the Kamioka mine in Japan 
and observes the electron anti- neutrinos emitted by ~ 16 nuclear power plants at distances 
of ~ 200 km. The results of KamLAND have provided convincing evidence for v^. disap- 
pearance and are in agreement with the so-called LMA-MSW solution of the solar neutrino 
problem (see, e.g.^ Ref. [24]). The neutrino oscillation analysis of current KamLAND data is 
dominated by statistical errors^ and it is and good approximation to gather various sources of 
systematical errors into an uncertainty on the overall number of events. However, in future, 
if more data are accumulated statistical errors will decrease and in principle one has to treat 
systematical errors more carefully. In this section we investigate the impact of uncertainties 
on the anti-neutrino flux for the determination of oscillation parameters in KamLAND. To 
this end we naively extrapolate the size of the data sample published in Ref. [6] to a total of 
five years data taking time by multiplying the event numbers by the factor 5 x 356/145.1, 
where 145.1 days is the exposure time of the reference sample. For further details on the 
KamLAND analysis see Refs. [24, 25] . 

In the case of KamLAND one has to generalize the expression for the number of events 
per positron energy bin from Eq. Q to account for the fact that several reactors (labeled by 
the index r) at different distances Lr contribute to the signal, and that neutrino oscillations 
occur: 

^« = ^Y. T^y.^ f dE,a{E,) MEu) ME,) Pee{Lr/E,) . (13) 
^ Ei J 

Here M is a, normalization constant, and frt are the power output and the element 
composition of the reactor r, respectively, and Pge is the oscillation probability depending 
on the neutrino mass squared difference Am? and the mixing angle 9. In our analysis we 
consider the following contributions to the covariance matrix of the event numbers Ni. 

V^"^ = NA, + N,N,al, + r,?- , (14) 

where the first term of the right hand side is the statistical error and the second term 
accounts for the over all normalization error a Act- The uncertainties on the anti-neutrino 
flux are included in V^^^^ which we take as 

Here 5Pr and 5fri are the errors on the power output and isotope composition of each 
reactor, and we assume typical values of {6Pr)/Pr = 0.02 and (Sfri)/ fri = 0.01. These errors 
are taken uncorrelated between the reactors. The last term in Eq. (fTH|) takes into account 
the uncertainty on the coefficients of the parameterization for the anti-neutrino fluxes as 
discussed in Sec. El 

In Figure El we show the 3a allowed regions for the oscillation parameters after 5 years of 
KamLAND data for various assumptions about the systematic errors. We observe from this 
figure that even after 5 years the KamLAND analysis is dominated by the statistical and 

■^The robustness of the KamLAND results with respect to statistical fluctuations has been extensively 
discussed in Ref. [25]. 
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Figure 3: (Color online) Scr allowed regions for sin'^ 29 and Am^ after 5 years of KamLAND data. The 
shaded region corresponds to statistical errors only, the regions delimited by the curves correspond to various 
assumptions about systematical errors. iTdet is a fully correlated error on the overall normalization. For the 
curves labeled "power, fuel, spectra" we include a 2% error on the power output of each reactor, 1% error 
on the fuel composition of each reactor, and the uncertainty on the anti-neutrino spectrum as described in 
the text. 



the overall normalization errors. The shaded region corresponds to statistical errors only, for 
the thin solid curve only the normalization error of a^et = 6.42% [6] is included. In Table 2 
of Ref. [6] various contributions to a^et are listed. If the uncertainties related to the flux are 
subtracted cxdet is reduced to 5.46% (see blue/dashed curve in Figure EI). For the red/solid 
and green/dash-dotted curves also the flux uncertainties according to Eq. (|TH|) are included. 
In both cases we find a very small effect on the oscillation parameters. 

To summarize, we find that even for 5 years of KamLAND data the determination of 
the oscillation parameters is dominated by statistical and overall normalization errors. The 
effect of flux uncertainties is rather small, and in particular it is not necessary to fully 
take into account flux shape errors. The inclusion of only the normalization errors for the 
total flux from each element (see Table ^ leads to nearly identical results as accounting 
for the full covariances of the coefficients a^- However, the proper treatment of the flux 
uncertainties reduces the overall normalization error. This will become more relevant if in 
future KamLAND analyses normalization errors like the uncertainty on the fiducial volume 
will be reduced. This will be relevant mainly for the measurement of the mixing angle, the 
determination of Am^ is hardly affected by any of the systematical errors. 

Finally, we note that already for the current KamLAND data sample [7] small differences 
in the allowed regions are visible due to the use of our neutrino fluxes, compared to the 
parameterization of Ref. [22]. For completeness, we mention that additional effects like the 
time evolution of the individual reactor powers or isotope compositions due to burn-up [26] 
may become relevant for future KamLAND analyses. The investigation of such effects is 
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Figure 4: (Color online) Difference between the positron spectra for no oscillations A^^^' and Nj^^^\ 
obtained by using the 3 parameter fit (Table O and 6 parameter fit (Table |3J to the anti- neutrino spectra, 
respectively. We assume a total number of events of 40000 and a typical isotope composition of ^35^ . 239p-^ . 
238u : 24ip^ ^ Q_57 . Q 3Q . Q Qg . Q Qg_ rpj^g ghadcd area corresponds to the la statistical error band, i.e., 

±y N'f^\ For 60 bins in positron energy and statistical errors only we find x^/dof — 139/60. 



beyond the scope of the present work. 



5 Application to future reactor experiments to measure ^13 

Let us now discuss the relevance of the flux uncertainties for reactor experiments planned 
to measure the small leptonic mixing angle ^13. It has been realized that the bound on this 
angle from previous experiments [4, 5] can be significantly improved if in addition to a far 
detector at a distance of order 2 km from the reactor a near detector at a few hundred meters 
is used. Due to the large number of events, the near detector provides accurate information 
on the reactor neutrino flux. Identical near and far detectors with normalization errors o"det 
below 1% will provide an accuracy on sin^ 26'i3 of order 0.01 (see, e.g., Refs. [9,27-30]). 

In Fig m we show the difference between the positron spectra predicted for no oscillations 
from the three and six parameter description of the neutrino spectra. The comparison 
with the statistical accuracy for a total number of events of 40000, which is typical for the 
far detector of these experiments [9], shows that for such experiments a precise model for 
the flux is needed. The small differences between the two parameterizations are clearly 
distinguishable by the statistical precision in the far detector: For 60 bins in positron energy 
we obtain xVdof = 139/60, where = T.iWt^^ - N^^^^f/Nf^^K Note that in the near 
detector is much worse, because of the larger number of events. 

In the following we adopt the six parameter model for the neutrino fluxes and investigate 
the impact of the errors of the coefficients on the sensitivity to the mixing angle. For 
simplicity we consider here only two-flavour neutrino oscillations characterized by the mixing 
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angle 9 and the neutrino mass-squared difference Am^. The number of events in a given 
positron energy bin i in the detector A {A = F) can be calculated by 

Nf = P^y^^ j dE, a{E,) MEu) Ri{E.) Pee{LA/E,) . (16) 
Similar to Ref. [27] we take into account the various systematical errors by writing 

Here the parameters a and describe the uncertainty on the reactor power and the detec- 
tor normalizations, respectively. The term QA^f accounts for the energy calibration (see 
Ref. [27] for details) and the last two terms in Eq. (fTTjl describe the uncertainty on the 
isotope fractions and the coefficients cm-, respectively. To test the oscillation parameters we 
consider the 

where the Tf^ depend on 6 and Am^, whereas Of^ = Nf-{9 = Am^ = 0). For each value of 9 
and Am^ Eq. (fTHjl has to be minimized with respect to the "pulls" a,bA, gA,Cey^ke- 

In Figure El we show the results of this analysis for a wide range of Am^. The values of 
Am^ relevant for the ^13 measurement are constrained by atmospheric neutrino data to the 
interval 1.4-10-3eV^ < Am^ < S.G-lQ-^eV^ at 3a (shaded region in Figure El). The regions of 
Am^ up to leV^ could be relevant for oscillations into hypothetical sterile neutrinos [31,32]. 
The solid curves in Figure El are calculated for perfectly known anti- neutrino flux, i.e., by 
fixing the coefficients ^ki in Eqs. (irTjl and (fTH|l to zero. For the dashed curves we use the 
errors obtained in the fit to the beta spectra as discussed in Sec. 01 For the dotted curves 
the anti-neutrino fiux coefficients are treated as free parameters in the fit, i.e., we drop the 
last term in Eq. p8|) when minimizing with respect to C,ke- In Figure El we consider two 
experimental configurations. One corresponds to an experiment with 6 ■ 10^ events in the 
far detector for no oscillations and near and far detector baselines of Lnd = 0.15 km and 
Lfd = 1.05km. This set-up is similar to the Double-Chooz proposal [30]. For the second 
configuration we have assumed a somewhat larger near detector baseline of Lnd = 0.7 km, a 
far detector baseline optimized for Am? ~ 2 ■ 10~^eV^ of Lnd = 1.7 km, and a rather high 
luminosity of 6 ■ 10^ events. 

Let us first discuss the region of Am^ relevant for the ^13 measurement. We observe from 
the figure that in the case of Lnd = 0.15 km the limit does hardly depend on the assumptions 
concerning the fiux uncertainty. Even the fiux free limit is not much worse than the limit for 
no error on the fiux, since at these short distances the near detector provides a very accurate 
determination of the anti- neutrino fiux. In contrast, the fiux uncertainty has some impact on 
the ^13 measurement if the near detector baseline is somewhat larger. In that case oscillations 
start to build up already between reactor and near detector and the uncertainties on the 
initial fiux become relevant (see also Figure 12 in Ref. [27]). In the region Am? > 5 • 10~^ eV^ 
the main information relevant for the limit is provided by the near detector. Hence the fiux 
uncertainties become even more relevant for both configurations in that region. We note 
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Figure 5: (Color online) The 90% CL limit on sin^ 29 from reactor neutrino experiments with near and 
far detectors as a function of Am^. The bound is shown for two experimental configurations as indicated 
in the figure. The solid curves correspond to no errors on the reactor neutrino flux, for the dashed curves 
the covariance matrix from the flt to the beta spectra is used, and for the dotted curves the coefficients for 
the anti-neutrino fiux are treated as free parameters in the fit. In all cases we have assumed CTo = 1% for 
the uncertainty of the reactor power, crdct = 0.6% for the detector normalization, cjca\ = 0.5% for the energy 
calibration, and aj — 1% for the error on the isotope fractions. The shaded region indicates the range of 
Am^ allowed at Scr from atmospheric neutrino data. 



that around Am^ ~ 2 (7) • 10~^ eV^ for the big (small) experiment the limit again becomes 
independent of the flux uncertainty. In that region rather fast oscillations occur at the near 
detector which still can be resolved by the detector, but cannot be mimicked by adjusting 
the coefficients of the fiux parameterization. Obviously, no limit can be obtained for the fiux 
free analysis in the averaging regime of very high Am^. 



6 Potential of an anti-neutrino detector close to a reactor 

In this section we investigate the potential of an anti-neutrino detector very close to a reactor, 
where "very close" is defined by the requirement that neutrino oscillations do not occur. This 
could be for example a near detector of the experiments considered in the previous section, 
if it is situated close enough to the reactor. In Subsec. 16.11 we investigate to what extent 
the uncertainty on the anti-neutrino fiux can be reduced by such a measurement, whereas in 
Subsec. 16.21 we consider the possibility to determine the isotope composition of the reactor 
core. 

To this aim we write the theoretical prediction for the number of events in bin i given in 
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Figure 6: (Color online) Improvement for the anti-neutrino flux errors for ^35^ function of the total 
number of events in a detector close to a reactor. We show S^™, where 5'"°* is defined in Eq. H21() and a 
runs over the 6 pulls associated with the flux- uncertainties shown in Fig. |2] We assume that the core contains 
97% ^^^U, and we take da — 1%, crdot — 0.6%, cJcai ~ 0.5%, and ct/ = 1%. The dotted hne corresponds 
roughly to the number of events expected in the near detector of the Double-Chooz experiment [30] . 

Eq. jnj) as 

T. = 0, + 5^p,$^„, (19) 

a 

where we drop the detector index A and we use the fact that Ni = Oi for no oscillations. 
The index a runs over all the pulls: Pa = {a, b, g, ^m) and the coefficients $^ can be read 
off from Eq. (fT7|) . With this notation Eq. ()18|) becomes 

i a ^ 

where 5pa is the error on the pull Pa-, which can be read off from Eq. pSj) . Departing from 
Eq. ()2()|1 it is straight forward to compute the improvement of the knowledge on a given 
parameter p^ due to the data Oi\ Because of the quadratic structure of Eq. the new 
covariance matrix S^^"^ of the Pa can be obtained by inverting 

Note that in general the Pa will be correlated after the measurement, i.e., 5'°'^™ will aquire 
non-diagonal entries from the second term in Eq. (PT|) . The final one sigma error on a 
parameter p^ is given by \/S^ . 

6.1 Improving our knowledge on the anti-neutrino flux 

From Figure one can see that the limit on the mixing angle is nearly the same for an 
analysis with completely free anti-neutrino fiux coefficients (dotted curves) and for the cur- 
rent errors on them (dashed curves), in the region Am^ < 10~^eV^, where oscillations can 
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be neglected in the near detector. This indicates that the near detector provides a rather 
precise determination of the flux on its own. 

This fact is quantified in Figure IHl where we show the improvement of the flux errors 
for 235u fj-om the near detector data with respect to the present errors obtained from the 
fit to the beta spectra measurements. We observe from the figure that for a number of 
events > 10^ the errors on the flux coefficients from the anti-neutrino measurement become 
comparable to the current errors. To reduce the errors by a factor two roughly 10^ events are 
needed. To avoid correlations with the flux coefficients from other isotopes, we assume that 
the core contains practically only ^^^U. However, for large number of events the errors on 
the coefficients Cki become strongly correlated. Especially the coefficients corresponding to 
the lines labeled "2" , "3" and "4" in Figure IHl are nearly fully correlated. This implies that 
a certain combination of these coefficients is severely constrained, and one should perform a 
diagonalization of the covariance matrix (similar as described in Section ^ to obtain again 
uncorrelated flux uncertainties in analogy to Figure |21 Note that the modes labeled "5" and 
"6" corresponding to relatively "fast oscillations" (compare Figure El) can be determined 
rather good by the anti- neutrino measurement. We find only modest correlations of these 
coefficients. 

In general also sizable fractions of ^^^Pu, ^^^Pu and ^ss^j present in the reactor. 

In this case all coefficients will become correlated. It might be possible to disentangle the 
contributions of the various isotopes taking into account precise information on the time 
evolution of the reactor composition. The improvement for a given isotope depends strongly 
on the relative amount of this isotope in the core. 

6.2 Determination of the isotope composition of a reactor 

In this subsection we investigate the possibility to determine the isotope composition of a 
reactor core by a nearby anti-neutrino detector. This could lead to applications of neutrino 
spectroscopy for reactor monitoring, either for improving the reliability of operation of power 
reactors or as a method to accomplish certain safeguard and non-proliferation objectives. In 
both cases the price tag of a moderately sized detector is small compared to the overall 
cost or benefit. Therefore the applicability of neutrino spectroscopy seems to depend only 
on its performance compared to existing technologies. For example the accuracy in the 
determination of the thermal power of civil power reactors as used for the production of 
electricity typically is in the range 0.6% — 1.5% [5]. The isotopic composition usually is 
not measured in situ but deduced from the time development of the reactor thermal power 
and the initial isotopic composition by using detailed reactor core simulation tools and is 
typically accurate at the percent level [11]. 

On the other hand, one requires for safeguard purposes to achieve a sensitivity which 
allows to quickly detect the diversion of weapons grade material at the level of one critical 
mass [10], e.g. for Plutonium this is approximately 10 kg. The average power reactor contains 
three tons of fissionable material of which roughly 40% have been converted to Plutonium 
by the end of the fuel life time, thus 10 kg Plutonium correspond to ~ 0.8% of the total 
content of Plutonium, which is beyond the accuracy levels for in situ monitoring today. The 
determination of the isotopic composition by the traditional methods furthermore relies on 
the assumption that the operator of the reactor collaborates. In the following we will show 
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Figure 7: (Color online) The Icr error on the fraction of ^^^U, ^^^Pu, and ^''^Pu as a function of the number 
of anti- neutrino events. The think curves correspond to the error on the sum of ^'^^Pu and ^^^Pu. The straight 
lines are calculated for perfectly known anti-ncutrino fluxes. For the sum of ^■^^Pu and ^''^Pu wc show also the 
result assuming the present anti-neutrino flux uncertainties ( "cTflux" ) , and uncertainties reduced by factors 
three ("crfiux/3") and ten ("crflux/10"). Also shown is the relative Icr error on the reactor power for crflux = 
and present flux errors. We assume an isotope composition ^35^ ; 239p^j . 238u . 24ip^ _ . . q j . g ^ 
and a detector normalization uncertainty (Jdet = 0.6%. 

that neutrino spectroscopy lias tlie potential to reach a sensitivity comparable to existing 
technologies and it does not require detailed information on the power history or the initial 
fuel composition, which is an important advantage especially in safeguards applications. In 
contrast to the existing literature on this topic [10-12] we use the full spectral information 
and therefore do not require an independent determination of the reactor power. In general, 
any safeguard regime based only on the total rate suffers from two problems: The first one is 
related to the availability of reliable information on the thermal power, whereas the second 
one is related to the fact that for most reactor types the diversion of core inventory is only 
possible during refueling, i. e. when the reactor is switched off and there is no neutrino flux. 
Thus in order to detect any diversion in this period an absolute measurement of the neutrino 
flux as well as of the thermal power is required. Moreover the composition of the new fuel 
has to be known exactly in order to predict the spectrum which is expected in case of no 
diversion. 

In the following we perform a fit where the isotope fractions are treated as free pa- 
rameters, subject to the condition ^ = 1. We impose no external information on the 
reactor power, i.e., no knowledge at all about the reactor is assumed. This means that we 
set l/{5paY = in Eq. (pT|l for a corresponding to cXa and cx/. The determination of the 
isotope fractions and the power is solely based on the differences between the anti-neutrino 
spectra emitted by the four isotopes. In Figure d we show the la accuracy obtained on the 
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Figure 8: (Color online) Contours of Ax^ — 1,4,9 in the plane of the ^'^^Pu and ^"^^Pu isotope fractions 
for 10^ anti-neutrino events. The colored regions correspond to perfectly known anti-neutrino flux shapes, 
whereas for the curves we assume uncertainties on the flux coefficients three times smaller than the present 
errors. We adopt the same isotope composition and add as in Figure [7| Reactor power and ^ssy fraction 
are treated as a free parameters. 

isotope fractions and the reactor power as a function of the anti-neutrino events. To give 
an example, for a detector with one ton fiducial mass at a distance of ten meters from a 
reactor with three GW thermal power roughly 10^ events are expected within 3 months of 
measurement time. 

First, one can see from Figure [7| that for > 10^ events a rather precise determination of 
the reactor power at the < 3% level is possible, given the current uncertainty on the anti- 
neutrino fluxes. For perfectly known fluxes the power accuracy is limited by the systematical 
uncertainty of the detector normalization. Second, for > 10^ events also the isotope fractions 
of ^^^U, ^^^Pu, and ^^^Pu can be determined at the percent level if no errors on the anti- 
neutrino fluxes are taken into account.^ The accuracy on the sum of the ^^^Pu and ^^^Pu 
fractions is clearly better than the one on the individual fractions. This is a consequence 
of the strong anti-correlation between the two Pu isotopes, which we illustrate in Figure |H1 
where contours in the ^^^Pu-^^^Pu plane are shown for 10® events. Note that for safeguard 
applications actually the sum of both Plutonium isotopes is the interesting quantity. 

From Figure m one can see that to determine the isotope composition a precise knowledge 
of the emitted fluxes is necessary. With present errors the la accuracy is limited to > 
10%. To reach a determination at the percent level the errors on the coefficients of the flux 
parameterization have to be reduced by a factor of three to ten. A factor three would be 
approximately achieved by the near detector of an experiment like Double-Chooz [30]. 

Let us note that in this analysis we do not take into account additional information such 
as the time evolution and reactor burn-up, or information from various traditional safeguard 

^We consider the fractions of ^^'"^U, ^^^Pu, and ^^^Pu as independent parameters, and determine the ^ssy 
fraction by the constraint = 1 • 
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methods. The main conclusion from the above results is that anti-neutrino spectroscopy 
may play an important role for reactor monitoring, especially since one expects significant 
synergies due to the combination with alternative technologies. 

7 Summary and conclusions 

In this work we have presented an accurate parameterization of the anti-neutrino flux pro- 
duced by the isotopes ^^^U, ^^^Pu, and ^^^Pu in nuclear reactors. Wc use a polynomial of 
order 5 and determine the coefficients by performing a fit to spectra inferred from experi- 
mentally measured beta spectra. Furthermore, the correlated errors on these coefficients are 
determined from the fit. 

Subsequently we investigate the impact of the flux uncertainties for the KamLAND ex- 
periment and ftitTirc reactor experiments to measure the mixing angle ^13. We show that 
flux shape uncertainties can be safely neglected in the KamLAND experiment, however the 
proper treatment of the errors associated to the anti-neutrino flux reduces somewhat the 
overall systematic error in KamLAND, which has some impact on the determination of the 
mixing angle. Future high precision reactor neutrino experiments with a far detector at dis- 
tances of order 2 km and a near detector at hundreds of meters are sensitive to the flne details 
of the reactor neutrino spectra. We find that a parameterization based on a polynomial of 
order two is not accurate enough to describe the anti-neutrino spectrum at the required level 
of precision. If the near detector is located at distances > 500 meters the flux uncertainties 
are relevant for the ^13 measurements. Moreover, in searches for sterile neutrinos at values 
of Am^ > 10^^ cV^ the main information is provided by the near detector, and hence the 
inclusion of anti-neutrino flux uncertainties is essential. 

Finally, we have investigated the potential of a detector very close to a reactor to improve 
on the uncertainties of the anti-neutrino fluxes, and to determine the isotopic composition 
in nuclear reactors through an anti-neutrino measurement. We flnd that without any ex- 
ternal knowledge on the reactor a three month exposure of a one ton detector allows the 
determination of the isotope fractions and the thermal reactor power at a few percent ac- 
curacy. This may open the possibility of an application for safeguard or non-proliferation 
objectives, which does not rely on information on the reactor thermal power or on the initial 
fuel composition, and hence neutrino spectroscopy can provide information complementary 
to traditional monitoring methods. To achieve this goal a reduction of the present errors 
on the anti-neutrino fluxes of about a factor of three is necessary, which naturally can be 
obtained from the data of the near detector of a Double-Chooz like experiment. 
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A Results of the fits to the anti-neutrino spectra 

In this appendix we give the best-fit coefficients for the polynomials used to parameterize 
the anti-neutrino flux of the isotopes ^^^U, ^^^Pu, and ^^^Pu according to Eq. In Table |21 
the coefficients for the polynomial of order 2 are given, wheres in Table 0] we display the 
coefficients, their errors and the correlation matrix for the polynomial of order 5. 





ai£ 


235 -Q 

239Pu 
24^PU 


0.904 -0.184 -0.0878 
1.162 -0.392 -0.0790 
0.852 -0.126 -0.1037 


238 -jj 


0.976 -0.162 -0.0790 



Table 2: Coefficients of tlie polynomial of order 2. For ^^^U, ^^^Pu, ^''^Pu the numbers are obtained from 
a fit to the data from Refs. [16, 17], for ^ssy reproduce the values given in Ref. [22]. 
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I = 235U 


correlation matrix p^^, 


k 




5aki> 


1 


2 


3 


4 


5 


6 


1 


3.519 


-10" 


7.26 


10-1 


1.000 


-0.996 


0.987 


-0.973 


0.956 


-0.938 


2 


-3.517 


- 10° 


8.81 


10-1 


-0.996 


1.000 


-0.997 


0.989 - 


-0.976 


0.962 


3 


1.595 


-10° 


4.06 


10-1 


0.987 


-0.997 


1.000 


-0.997 


0.990 


-0.980 


4 


-4.171 ■ 


10-1 


8.90 


10-2 


-0.973 


0.989 


-0.997 


1.000 - 


-0.998 


0.992 


5 


5.004 ■ 


10-2 


9.34 


10-3 


0.956 


-0.976 


0.990 


-0.998 


1.000 


-0.998 


6 


-2.303- 


10-3 


3.77 


10-^ 


-0.938 


0.962 


-0.980 


0.992 - 


-0.998 


1.000 






= 239Pu 




correlation matrix p^^, 


k 




5aki> 


1 


2 


3 


4 


5 


6 


1 


2.560 


-10" 


4.01 


10-1 


1.000 


-0.993 


0.977 


-0.954 


0.928 


-0.899 


2 


-2.654 


-10° 


5.58 


10-1 


-0.993 


1.000 


-0.995 


0.982 - 


-0.962 


0.938 


3 


1.256 


■10« 


2.91 


10-1 


0.977 


-0.995 


1.000 


-0.996 


0.984 


-0.967 


4 


-3.617- 


10-1 


7.17 


10-2 


-0.954 


0.982 


-0.996 


1.000 - 


-0.996 


0.986 


5 


4.547- 


10-2 


8.37 


10-3 


0.928 


-0.962 


0.984 


-0.996 


1.000 


-0.997 


6 


-2.143 - 


10-3 


3.73 


10-^ 


-0.899 


0.938 


-0.967 


0.986 - 


-0.997 


1.000 






= 241pu 




correlation matrix p^^, 


k 






1 


2 


3 


4 


5 


6 


1 


1.487 


-10° 


3.23 


10-1 


1.000 


-0.991 


0.974 


-0.950 


0.923 


-0.893 


2 


-1.038 


-10° 


4.31 


10-1 


-0.991 


1.000 


-0.994 


0.980 - 


-0.960 


0.936 


3 


4.130- 


10-1 


2.15 


10-1 


0.974 


-0.994 


1.000 


-0.995 


0.984 


-0.966 


4 


-1.423 - 


10-1 


5.02 


10-2 


-0.950 


0.980 


-0.995 


1.000 - 


-0.996 


0.986 


5 


1.866- 


10-2 


5.54 


10-3 


0.923 


-0.960 


0.984 


-0.996 


1.000 


-0.997 


6 


-9.229 - 


10-^ 


2.33 


10-* 


-0.893 


0.936 


-0.966 


0.986 - 


-0.997 


1.000 



Table 3: Coefficients ake of tlie polynomial of order 5 for the anti-neutrino flux from elements £ = ^^^\J, 
^^^Pu, and ^'^^Pu. In the column Sate the la errors on ake are given. Furthermore the correlation matrix of 
the errors is shown. 
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